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Abstract 

Dynamical linked cluster expansions are linked cluster expansions with 
hopping parameter terms endowed with their own dynamics. This amounts 
to a generalization from 2-point to point-link-point interactions. We develop 
an associated graph theory with a generalized notion of connectivity and de- 
scribe an algorithmic generation of the new multiple-line graphs . We indicate 
physical applications to spin glasses, partially annealed neural networks and 
SU(N) gauge Higgs systems. In particular the new expansion technique pro- 
vides the possibility of avoiding the replica-trick in spin glasses. We consider 
variational estimates for the SU(2) Higgs model of the electroweak phase tran- 
sition. The results for the transition line, obtained by dynamical linked cluster 
expansions, agree quite well with corresponding high precision Monte Carlo 
results. 
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1 Introduction 



Linked cluster expansions (LCEs) have a long tradition in statistical physics. Orig- 
inally applied to classical fluids, later to magnetic systems ([0,@,@ and references 
therein), they were generalized to applications in particle physics in the eighties [|J. 
There they have been used to study the continuum limit of a lattice $ 4 field theory 
in 4 dimensions at zero temperature. In [| they were further generalized to field 
theories at finite temperature, simultaneously the highest order in the expansion 
parameter was increased to 18. Usually the analytic expansions are obtained as 
graphical expansions. Because of the progress in computer facilities and the devel- 
opment of efficient algorithms for generating the graphs, it is nowadays possible to 
handle of the order of billions of graphs. The whole range from high temperatures 
down to the critical region becomes available, and thermodynamic quantities like 
critical indices and critical temperatures are determined with high precision (the 
precision is comparable or even better than in corresponding high quality Monte 
Carlo results) 0-0. An extension of LCEs to a finite volume in combination with 
a high order in the expansion parameter turned out to be a particularly powerful 
tool for investigating the phase structure of systems with first and second order 
transitions by means of a finite size scaling analysis [[10 . 



Linked cluster expansions are series expansions of the free energy and connected 
correlation functions about an ultralocal, decoupled theory in terms of a hopping 
parameter K. The corresponding graphical representation is a sum in terms of 
connected graphs. The value of K parametrizes the strength of interactions between 
fields at different lattice sites. Usually they are chosen as nearest neighbours, but also 
more general, less local couplings lead to convergent expansions (|Tl[ and references 
therein) under appropriate conditions on the decay of the interactions. In contrast 
to the ultralocal terms of a generic interaction we will sometimes refer to hopping 
terms as non-ultralocal. 

In this paper we develop dynamical linked cluster expansions (DLCEs). These 
are linked cluster expansions with hopping parameter terms that are endowed with 
their own dynamics. Such systems are realized in spin glasses or partially annealed 



neural networks with (fast) spins and (slow) interactions [|T^]-[|T4|. They also occur 
in variational estimates for SU(N)-gauge-Higgs systems as we will show later in this 
paper. Like LCEs they are expected to converge for a large class of interactions. 

Formally DLCEs amount to a generalization of an expansion scheme from 2-point 
to point-link-point-interactions. These are interactions between fields associated 
with two points and with one pair of points called link. The points and links are not 
necessarily embedded on a lattice, and the links need not be restricted to nearest 
neighbours. We develop a new multiple-line graph theory in which a generalized 
notion of connectivity plays a central role. Standard notions of equivalence classes of 
graphs like 1-line irreducible and 1- vertex irreducible graphs have to be generalized, 
and new notions like 1-multiple-line irreducible graphs must be defined in order to 
give a systematic classification. We describe algorithms for generating these classes. 
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This is essential, since the number of graphs rapidly increases with the order in the 
expansion parameter. It typically exceeds the number of LCE-graphs by more than 
an order of magnitude. 

The paper is organized as follows. In section 2 we specify the models that admit 
a DLCE. We introduce multiple-line graphs and explain the idea behind the abstract 
notions of multiple-line graph theory. Detailed definitions of multiple-line graphs , 
related notions and the computation of weights are given in section 3. Section 4 
treats the issue of renormalization in the sense of suitable resummations of graphs. 
Algorithms for their generation are described in section 5. Applications to spin 
glasses and (partially annealed) neural networks are indicated in section 6. There it 
is of particular interest that DLCEs allow for the possibility of avoiding the replica 
trick. In section 7 we present results for the transition line of an SU(2) Higgs 
model. These are obtained by DLCEs applied to gap equations that follow from 
convexity estimates of the free energy density. The results are in good agreement 



with corresponding high precision Monte Carlo results ||15|| . Section 8 contains a 
summary and outlook. 

The reader who is more interested in our actual applications of DLCEs may skip 
sections 3-5 in a first reading. These sections are of interest from a systematic point 
of view and essential if he is interested in applying DLCEs. 



2 A Short Primer to DLCEs 

In this section we first specify the class of models for which we develop dynamical 
linked cluster expansions. Next we illustrate some basic notions of multiple-line 
graph theory, in particular the need for a new notion of connectivity. 

By A we denote a finite or infinite set of points. One of its realizations is a 
hypercubic lattice in D dimensions, infinite or finite in some directions with the 
topology of a torus. Ai denotes the set of unordered pairs (x,y) of sites x,y G A , 
x ^ y, also called unoriented links, and Ai a subset of A^ 

We consider physical systems with a partition function of the generic form 

Z(H,I,v) = expW(H,I,v) 

= N ( V(f)VUexp (-S(<p, U, v)) exp ( H(x)<p{x) + ^ I(l)U(l)), (1) 

J A . -T- 



with measures 



and action 



V(f) = Y[ d(f)(x) , VU = Y[ dU{l) (2) 

zeA ieAi 



S(<j>, U,v) = + <^v)4>{x)U{x,y)4>{y\ (3) 
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with non-ultralocal couplings 



v(x,y) = v(y,x) ^ only for (x, y) E A 1 , 

in particular v(x, x) = 0. (4) 

For later convenience the normalization via M is chosen such that W[0, 0, 0] = 0. 

The field <p(x) is associated with the sites x G Ao and the field U(l) lives on the 
links I G Ai, and we write U(x,y) = U(l) for I = (x,y). For definiteness and for 
simplicity of the notation here we assume <f>(x) G R and U(l) G R. In our actual 
applications is a 4-component scalar field and U an SU(2)-valued gauge field, or 
and U both are Ising spins, or the 0s are the (fast) Ising spins and the Us G R 
the (slow) interactions. The action is split into two ultralocal parts, S° depending 
on fields on single sites, and S 1 depending on fields on single links I G A x . For 
simplicity we choose S 1 as the same function for all links I G Ai. We may identify 
Ai with the support of v, 

Ai = {l = (x,y)\v(x,y)^0}. (5) 

The support of v(x,y) need not be restricted to nearest neighbours, also the precise 
form of S° and S 1 does not matter for the generic description of DLCEs, S° and 
S 1 can be any polynomials in and U, respectively. The only restriction is the 
existence of the partition function. In one of our applications S° will be a 4 -type 
theory with O(N) symmetry, S l consists of a term linear or quadratic in U. 

Note that the interaction term v(x,y) 0(x) U(x,y) <p(y) contains a point-link- 
point-interaction and generalizes the 2-point-interactions v(x,y) <fi(x) 4>(y) of usual 
hopping parameter expansions. The effective coupling of the fields has its own 
dynamics governed by 5' 1 (?7), the reason why we have called our new expansion 
scheme dynamical LCE. 

Dynamical linked cluster expansions are induced from a Taylor expansion of 
W(H, I, v) = In Z(H, I, v) about v = 0, the limit of a completely decoupled system. 
We want to express the series for W in terms of connected graphs. Let us consider 
the generating equation 

dW/dv(xy) = 1/2 < <j>(x)U(x,y)<j>(y) > 

= 1/2 (W H y)H{y) + W H{x)H ( y) W I{ 

x ,y) 

+ W H{x) W H{y) W I{x , y) ). (6) 

Here < ■ > denotes the normalized expectation value w.r.t. the partition function 
of Eq. (H). Subscripts H(x) and I(x,y) = I(y,x) = 1(1) denote the derivatives of 
W w.r.t. H(x) and I(x,y), respectively. 

Next we would like to represent the right hand side of Eq. (|6|) in terms of con- 
nected graphs. Once we have such a representation for the first derivative of W 
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w.r.t. v, grapical expansions for the higher derivatives can be traced back to the 
first one. 

For each W in Eq. (|6|) we draw a shaded bubble, for each derivative w.r.t. H 
a solid line, called a 0-line, with endpoint vertex x, and for each derivative w.r.t. 
I a dashed line, called a [/-line, with link label I = (x,y). The main graphical 
constituents are shown in Fig. 1. Two 0-lines with endpoints x and y are then 
joined by means of a dashed [/-line with label I, if the link I has x and y as its 
endpoints, i.e. / = (x,y). According to these rules Eq. (0), multiplied by v(x,y) and 
summed over x and y, is represented by Fig. 2. Note that, because of the Taylor 
operation, each solid line from x to y carries a factor v(x,y). 

Since the actual need for a new type of connectivity is not quite obvious from 
Fig. 2, because Eq. (^j does not contain higher than first order derivatives w.r.t. /, 
let us consider a term 



ii 




In 

Figure 1: Graphical representation of the derivatives of W(H,I,v). (a) n-point 
function d n W/dH{x 1 ) ■ ■ ■ dH(x n ), (b) n-link function d n W/dI(h) ■ ■ ■ dl(l n ). 




Figure 2: Generating equation of the graphical expansion of DLCEs. The solid line 
in each graph carries a propagator v(x, y). A dashed [/-line with label I intersects a 
solid line with endpoints x and y if I = (x, y). 

W H[x) W H (y) W H(r) W H (s) W I{x>y)I{ r,s) ( 7) 
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occurring in the second derivative of W w.r.t. v(x,y), v(r,s). According to the 
above rules this term would be represented as shown in Fig. 3a. While the 2 
vertices in the last term of Fig. 2 are connected in the usual sense via a common 
(solid) line (the dashed line with an attached bubble could be omitted in this case), 
the graph in Fig. 3a would be disconnected in the old sense, since neither x nor y 
are line-connected with r and s, but -as a new feature of DLCE graphs- the lines 
from x to y and from r to s are connected via the dashed lines emerging from a 
common bubble shown in the middle of the graph. As we see from Fig. 3a, we 
need an additional notion of connectivity referring to the possibility of multiple-line 
connectivity. While the analytic expression is fixed, it is a matter of convenience 
to further simplify the graphical notation of Fig. 3a at v — 0. Two possibilities are 
shown in Fig. 3b and Fig. 3c. To Fig. 3b we later refer in the formal definition 
of the new type of multiple-line connectivity. In the familiar standard notion of 
connectivity two vertices of a graph are connected via lines. The vertices are line- 
connected. Already there, in a dual language, one could call two lines connected 
via vertices. The second formulation is just appropriate for our need to define when 
two lines are connected. The corresponding vertices mediating the connectivity of 
lines are visualized by tubes, in Fig. 3b we have just one of them. The tubes should 
be distinguished from the former type of vertices represented as full dots which are 
connected via bare ^-lines. In Fig. 3c we show a simplified representation of Fig. 3b 
that we actually use in graphical expansions. 




Figure 3: Representation of WH( x )WH(y)WH( r )WH(s)Wi(x,y)i(r,s)- ( a ) according to 
the rules of Fig. 1 and 2, (b) same as (a), but at v — and simplified for a formal 
definition of multiple-line connectivity, cf. section 3, (c) same as (b), but for use in 
the actual graphical representations. 

The derivative terms have to be evaluated at v = 0. For v = we have a 
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decomposition of W according to 



W(H, I,v = 0)=J2 W°(H(x)) + W\l(l)) (8) 



with 



and 



expW(ff) ee Z\H) = j - J p '.UK (9) 

J_ oo c/0exp(-5°(0)) 

expW^l) = Z 1 (/) = ^_ yy y ! , (io) 

In Eq.s (|],|10]) we have omitted any single site or single link dependence, because 
we assume that 5° and S l are the same for all x G A and I G A 1; respectively. 
Therefore, at v — 0, the only non-vanishing derivatives of W are 



a n iy°(jj-(o;i )) 



H / H(a;i)H(a;2)...fl'(a!n)| w= o _ f>uY„ \ n Ox x ,x 2 ,.~xn l iJ -J 




and 

^/(h)/^).../^)!^ - — di(h) m — ° h ' h >- l ™> { > 

but mixed derivatives w.r.t. H and I vanish. As anticipated in Fig.s 3b and 3c, for 
v = we replace the dashed bubbles and graphically distinguish between bubbles 
with Klines and bubbles with [/-lines. We define 

for a connected n-point vertex with n> 1 bare 0-lines emerging from it and 

" = m " = { ~ai^),. (14) 

for a connected z/-line consisting of v bare lines. If v — 1, we often omit the dashed 
line. If the bare lines of a z/-line are internal 0- lines, they get vertices attached to 
their endpoints, if they are external [/-lines, no vertices will be attached. 

Let V denote the lattice volume in D dimensions. The Taylor expansion of W 
about v = to second order in v then reads 

W(H,I,v) = W(H, I,v = Q) 

+ ^2 v(x,y) -W H(x) W H(y) W I{x , y) 

x,y€Ao 
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+ 2 ^(^M r ,s) 

x,j/,r,sGAo 



+ 2 W / ff(z) J H'(r)W / H(y)H(s) W^(j;,!/)Wj( r)8 ) 

+ 4 W^Hfy) W^ff (s) W^ff^)//^) Wl{x,y)I{r,s) 

+ 2 W / H(r) J H"(x)W / H(y) J H"(s)W / 7(r,y)/(a;,s) 

+ Wh(x) Wh( v ) W H (r) Wh(s) ^J(i,j);(r,s) 

+ 0(v 3 ), 



v=0 



(15) 



where we have used that v (x, x) = 0. For each W in the products of Ws we now 



insert Eq.s ([[]]), (0). 

If we choose v in a standard way as next-neighbour couplings 



D-l 



v(x, y) = 2K y^((5 x+A , y + 5 : 

Ai=0 



(16) 



with p, denoting the unit vector in /^-direction, Eq. fllBD becomes in a graphical 
representation at H — I — 



W{0,0,v) 
V 



= (2K) - (2D) - 
+ (2iT) 2 {i (2D) 2 



+ 2 (2^) 2 



+ - 2(2L>) 




+ ^ (2D) 




(17) 



+ 0(K 3 ). 

For clarity, here we have written explicitly the topological symmetry factors and 
the lattice embedding numbers. (Usually graphs represent their full weights includ- 
ing these factors.) Note that the first two graphs of the second order contribution 
also occur in a usual LCE with frozen [/-dynamics, the next two differ by an addi- 
tional dashed 2-line and the last one becomes even disconnected without the dashed 
line. 

As usual, graphical expansions for correlation functions, in particular suscepti- 
bilities, are generated from W(H, J, v ) by taking derivatives w.r.t. the external fields 



S 



H and /. Graphically this amounts to attaching external 0-lines and ?7-lines with 



'1 endpoint) attached to vertices, e.g. 




(no endpoint) attached to i^-lines, e.g. 



In passing we remark that the conventional LCE is included as a special case of the 
DLCE, if the [/-dynamics is "frozen" to some value Uo 7^ 0, so that 



W\I) 
dW\l) 
dl 

dl n 



U , 



IUa 



(19) 



for all n > 1, 



i.e., no n-lines do occur with n > 1. In this case it becomes redundant to attach 
dashed lines to bare lines. As mentioned above, in an LCE only the first three 
contributions would be left in Eq. (|TTD. 



3 Graphical expansion 
3.1 Multiple-line graph theory 

The definition of a multiple-line graph as it will be given here is adapted to the 
computation of susceptibilities, where the sum is taken over all possible locations of 
the fields. The definition easily generalizes to correlation functions. 

For details of the standard defmiton of graphs in the framework of linked cluster 
expansions and related notions we refer e.g. to 0, Here, for convenience, we briefly 
recall the very definition of a graph to point out the new properties of multiple-line 
graphs as defined below in this section. 

A (standard LCE) graph or diagram is a structure 

f = (£ T ,Br,Er,$r), (20) 

where £r and Br ^ are disjoint sets of internal lines and vertices of T, respectively. 
Er is a map 

E r :B r - {0,1,2,...}, 

v -> E T {v) (21) 

that assigns to every vertex v the number of external lines Er(v) attached to it. 
Finally, $r is the incidence relation that assigns internal lines to their two endpoints. 
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A multiple-line graph or multiple-line diagram is a structure 

r = (Cr,Mr,B r ,4*\4 U) ,*r,*r). (22) 
Cr, M-t and Br are three mutually disjoint sets, 

Cr = set of bare internal lines of T, (23) 

Mr = set of multiple lines of T, (24) 

Br = set of vertices of T. (25) 

E r ' is a map 

4*>:Br - {0,1,2,...}, 

u -> 4%) (26) 

that assigns to every vertex v the number of bare external 0-lines E r ^\v) attached 
to v. Every such 0-line represents a field 0. The number of external 0-lines of T is 
denoted by = £„ eBr E^\v). Similarly, E T U) is a map 

EP:Mr - {0,1,2,...}, 

m -> ^fV) (27) 

that assigns to every multiple line m the number of external [/-lines E r U \m) at- 
tached to m. Every such [/-line represents a field U associated with a lattice link. 
The number of external [/-lines of T is given by E T U ^ = J2 meMr E T u \m). 

Furthermore, $r and \l/r are incidence relations that assign bare internal lines 
to their endpoint vertices and to their multiple lines, respectively. We treat lines as 
unoriented. The generalization to oriented lines is easily done. More precisely, let 
(Br x Br) be the set of unordered pairs of vertices (v,w) with v,w G Br, v ^ w. 
(The bar implies unordered pairs, the prime the exclusions of (v,v), v G Br-) As in 
the standard linked cluster expansion, self-lines are excluded. Every bare internal 
line is then mapped onto its pair of endpoints via 



$r : C r -> (Br x Br) ■ (28) 

We say that v and w are the endpoint vertices of I G Cr if $r(0 = (v,w). If there 
is such an Z G Cr, v and w are called neighbours. Similarly, \l/ r is a map 

I -> * r (Z) (29) 

that maps every bare internal line to a multiple line. A multiple line m G A^r 
is composed of bare internal lines I G Cr which belong to m in the sense that 
\l/r(0 = Tn. Ix r (m) is the total number of bare internal lines belonging to m. With 
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v = lM T ( m ) + -^t ( m )> m i s called a z/-line. We always require that v > 1. On the 
other hand, every bare internal line belongs to one and only one multiple line. For 
simplicity we often identify a 1-line with the only one bare line that belongs to it. 

Next we introduce some further notions that will be used later. External vertices 
are vertices having external 0-lines attached, 

B r , ext = {veB r \El?\v)^0}, (30) 

whereas internal vertices do not, Br^ n t — Br\ Br^xt- Similarly, external multiple 
lines have external [/-lines attached, 

Mr,ext = {meM r |4 f/) H^}» (31) 

and the complement in Air are the internal multiple lines, Air, int = Air \ Mr,ext- 
For every pair of vertices v,w G Br, v ^ w, let <£> (v,w) be the set of lines with 
endpoint vertices v and w, and |$ (v,w)\ the number of these lines. Thus $ (v,w) 
is the set of lines v and w have in common. With E^\v) denoting the number of 
external 0-lines attached to v G Br, 

t Br {v) = £ \*\v,w)\ + 4%) (32) 

«>e£>r 

is the total number of bare lines attached to v. 

Some topological notions and global properties of graphs will be of major interest 
in the following. A central notion is the connectivity of a multiple-line graph . Recall 
that we want to consider the DLCE expansion of the free energy and of truncated 
correlation functions as an expansion in connected graphs. As indicated in section 
2, the main generalization compared to the common notion of connectivity of a 
graph which is required here is that an additional type of connectivity is provided 
by multiple-lines. To define the connectivity of a multiple-line graph T, T first is 
mapped to a (standard) LCE graph T to which the standard notion of connectivity 
applies. There are various equivalent ways to define such a map. We choose the 
following one. 

• For every multiple-line m G Air define a new vertex w(m). Let Br = 
{w(m)\m G Air} an d define B = Br U Br as the union of the vertices of 
r and the new set of vertices originating from the multiple-lines. 

• For every bare internal line / G Cr define two new internal lines /i,/ 2 and 
incidence relations 

£(Zx) = (v^w^ril))), 

$(/ 2 ) = foMMQ)), (33) 

where v\ and v 2 are the two endpoint vertices of /. The set of all lines h,h, 
for all / G Cr, is denoted by C 
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• Define the external incidence relations 



E : B - {0,1,2,...}, 

E(v) = E&\v), forveB r , (34) 
E{v) = E^\m), for v = w(m) G £>r- 

Now, T is defined by 

f = (C,B,E,^). (35) 

Having defined the standard LCE graph T for any multiple-line graph T, we call 
T multiple-line connected or just connected if T is connected (in the usual sense). In 
Fig. 4 we have given two examples for a connected (upper graph) and a disconnected 
(lower graph) multiple-line graph . 





Figure 4: Example of multiple-line connectivity. The upper multiple-line graph Ti 
is connected because the graph Ti is connected in the conventional sense. The lower 
multiple-line graph T 2 is disconnected because T 2 is so. 

The next important notion is the topological equivalence of two multiple-line 
graphs . Two multiple-line graphs 

r, = (CuMuB^E^ ,EP,^,%), i = l,2 (36) 

are called (topologically) equivalent if there are three invertible maps 

0i:#i - B 2 , 

02 : A -> A>, (37) 



such that 



$ 2 o 2 = X O $ 1; 

02 = 03°^1, (38) 
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and 



E^o^ = Ef\ (39) 
Here o means decomposition of maps, and 



4> 1 :B 1 xB 1 -> B 2 x B 2 

Mv,w) = (Mv)Ai(w)). (40) 

A symmetry of a multiple-line graph T = (£, A4, B : E^\ E^ u \ $, \1/) is a triple 
of maps </>i : B — > £>, 02 : £ — ► £ and 03 : .M — > .M such that 

$ o 2 = 1 o $, 

\pO02 = 03 O \1>, (41) 

and 

£?W o 0! = 

EWofa = E^\ (42) 

The number of these maps is called the symmetry number of V. 

We denote by Gei,e 2 (L) the set of equivalence classes of connected multiple- 
line graphs with L bare internal lines, E\ external 0-lines and E 2 external [/-lines. 
Furthermore we set 

G El ,E 2 := |J Ge x ,e 2 {L). (43) 

L>0 

A multiple line graph V does not need to have a vertex. If Br = 0, we have 
Cr = as well. If in addition T is connected, Air consists of only one element, 
with all external [/-lines attached to it. (We anticipate that T is 1-multiple-line 
irreducible (1MLI) by definition. For the definition of 1MLI cf. section 4 below.) 
The only graph of Go,e(L = 0) is given by 



E. (44) 



It represents the leading term of the susceptibility 



X0,E VD 

h,...,l E eAi 



l - <U(h)---U(l E )> c (45) 



and is given by d W (I)/dI | . The index c in (pE5|) stands for truncated (con- 
nected) correlation. 
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By removal of a z/-line m G Air we mean that m is dropped together with all bare 
internal lines and all external [/-lines that belong to m. This notion is explained in 
Fig. 5a. (It is used in section 4 for 1-lines to define 1-particle irreducible (1PI) and 
1-line irreducible (1LI) multiple-line graphs .) 

On the other hand, by decomposition of a u-line m G A4r we mean that m is 
dropped together with the external [/-lines of m, but all bare internal lines that 
belong to m are kept in the graph, being identified now with 1-lines. This notion 
will be used below to define 1MLI and renormalized multiple-line moments. It is 
illustrated in Fig. 5b. 

Similarly, decomposition of a vertex v G Br means to remove the vertex v and 
to attach the free end of every line that entered v before to a new vertex, a sepa- 
rate one for each line. This notion is used to define 1- vertex-irreducible (1VI) and 
renormalized vertex moments for multiple-line graphs . For an example see Fig. 5c. 




i 



Figure 5: Removal (a) and decomposition (b) of a 5-line. Decomposition (c) of a 
vertex. 



3.2 Susceptibilities and weights 

In the last section we have defined multiple-line graphs and the notions of connec- 
tivity and equivalence of such graphs. The definition is chosen in such a way that 
the series expansions of the free energy and of truncated correlation functions are 
obtained as a sum over equivalence classes of connected multiple-line graphs . The 
number L of bare internal lines of a multiple-line graph T counts the order in the 
expansion parameter v(x,y) to which T contributes. If v(x,y) is of the form 

v{x,y) = 2K J2 S v,*> ( 46 ) 
zeJ\f(x) 
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with J\f(x) any finite x-dependent set of lattice sites, the contribution of T is a 
multiple of {2K) L . Often used special cases are the nearest neighbour interactions 

D-l 

v(x, y) = (8x, y +ii + S x , y - n ) (47) 

and the uniform interaction 

v(x,y) = 2K (1-4,,), (48) 
which is used in models of spin glasses and partially annealed neural networks. 
Susceptibilities of the and U fields will be represented as 

Xe u e 2 = ^ E E <<P{x 1 )---<p(x El )U{l 1 )---U{l E2 )> c 

x 1 ,...,x El £A l 1; ...,l E2 eAi 
J_ x ^ d E ^W(H,I, V ) 

= ™ xl ,XeA 0h ,± e - M OHM • -m^dKk) • ■■dI(l E2 ) 

= E( 2 ^) L E w ^ ( 49 ) 

l>o reg El ,E 2 (L) 

with lattice volume V and dimension D. Similar representations hold for higher 
moments fi. 

The weight w(T) of a multiple-line graph T G Qe^e-JJ-^) * s gi ven as the product 
of the following factors 

• for every vertex v G Br a factor 

V - = \—dH^) H ^ (50) 
where n = tj3 r (v) is the total number of bare lines attached to v. 

• for every multiple line m G M.r a factor 

,1c . /'PWV) 



H=I=0 



l-flzHw (51) 
where z/ = lM r ( m ) + ^r\' m )i that is m is a z/-line, 

• a factor 1/Sr, where Sp is the topological symmetry number of T, 

• a factor counting the permutation symmetry of external 0-lines, 
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(52) 



a factor counting the permutation symmetry of external [/-lines, 

E m < 



rwV'M! ' (53) 

• the lattice embedding number of T, which is the number of ways T can be 
embedded on a lattice of given geometry, e. g. on a hypercubic lattice. To this 
end, the vertices of T (if any) are placed onto lattice sites. One arbitrary vertex 
is placed at a fixed lattice site, in order to account for the volume factor 1 /V 
in P9|). A priori there is no exclusion principle. This means that any number 
of vertices can be placed at the same lattice site. (This is sometimes called free 
embedding.) Two restrictions apply to the embeddings. The first constraint 
results from the fact that a bare internal line represents a hopping propagator 
v(x, y), with lattice sites x and y at which the two endpoint vertices of the line 
are placed at. A reasonable computation of the embedding number takes into 
account the particular form of v(x,y) from the very beginning. The second 
constraint is that bare lines of the same multiple-line have to be mapped on 
the same pair of sites. 

For example, if v(x,y) is the nearest neighbour interaction P7|), two vertices 
which have at least one line in common are to be placed at nearest neighbour 
lattice sites. On the other hand, a propagator v(x,y) of the form (|48| ) implies 
a rather weak constraint in that x and y must be different, but otherwise can 
be freely placed over the lattice. 

We remark that in case of a non-trivial internal symmetry (such as considered 
in section 7) the expressions of Eq.s (|49|)-(|5lD must be modified appropriately. In 



particular, the weight (51) of a multiple- line does no longer take such a simple form. 
Eq. ( |1 25|) below is an example for the case of a hopping term originating in an SU(2) 
Higgs model. 



4 Renormalization 

Truncated correlation functions, susceptibilities and other moments are obtained as 
sums over multiple-line graphs that are connected. Their number rapidly grows with 
increasing order, that is with increasing number of bare internal lines. The procedure 
of "renormalization" means that the connected moments are represented in terms 
of reduced ones. The reduced moments are obtained by summation over multiple- 
line graph classes which are more restricted than just by their property of being 
connected. Of course the number of graphs of such classes is smaller. Only the 
most restricted multiple-line graph classes must be constructed. The subsequent 
steps towards the moment computation are most conveniently done by operating 
analytically with the reduced moments. In particular, it is no longer necessary 
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to generate all connected and the corresponding intermediate multiple-line graph 
classes. 

A connected multiple-line graph L is called 1-particle irreducible (1PI) if it sat- 
isfies the following condition. Remove an arbitrary 1-line of L. There is at most one 
connected component left that has external lines attached. (This notion is the same 
as the one used in the context of Feynman graphs.) On the other hand, if in addition 
the remaining graph is still connected, then L is called 1-line irreducible (1LI). In 
many cases it is sufficient to use only the second notion. It is for instance sufficient 
that all vertices are constrained to have only an even number of lines attached, or 
more generally, if graphs and subgraphs with one external line are forbidden. For 
notational simplicity we assume in the following that this is the case and henceforth 
refer only to the notion 1LI []. The generalization to the case in which 1LI and 1PI 
graphs must be distinguished goes along the same lines as for LCEs, which was 
discussed in 

By GeYe 2 (.-L) we denote the subset of multiple-line graphs L 6 Qe 1 ,e 2 (L) that 
are ILL lLI-susceptibilities are defined as series in the hopping parameter similarly 
as in P9| ) by restricting the summation to 1LI graphs, 

Susceptibilities are easily obtained in a closed form in terms of lLI-susceptibilities 
X 1LI ■ It can be shown that the x 1LI s can De obtained by an appropriate Legendre 
transform. For instance 

xl L i 

(i - v{o)x2% r 

where v(k) is the Fourier transform of the hopping propagator v(x, y), 

/T rlDU 
_ e-*<*-*) «(*). (56) 

In LCEs the second important resummation comes from so called vertex renor- 
malizations. This means partial resummation of graphs with specific properties such 
as having one external vertex only. These sums then are considered as " renormalized 
vertices" replacing the vertices of graphs with complementary properties. The pro- 
cedure naturally leads to the notion of 1-vertex irreducibility (1VI) and renormalized 
moments. 

In DLCE we follow this procedure. The very definition of 1VI has to be modified 
slightly for multiple-line graphs because of the enhanced connectivity properties due 



L In ref. j|. || the term 1PI was used instead. 
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to multiple-lines. In addition, as a natural generalization, we supplement vertex 
renormalization by multiple-line renormalization. 

A multiple-line graph T is called 1- vertex irreducible (1VI) if it satisfies the 
following condition. Decompose an arbitrary vertex v G Br- Every connected 
component of the remaining graph has then at least one external line attached. It 
can be a 0-line or a [/-line. We write 

e%k(L) = {r g g^ E2 (L) | r is ivi} (57) 

for the set of equivalence classes of graphs that are both 1LI and IVI, with E\ 
external 0- lines, E 2 external [/-lines and L bare internal lines. 

The renormalized vertex moment graphs are 1LI graphs that have precisely one 
external vertex and no external multiple line, 

Q k (L) = {r G Gl${L) | there is v G B T with E^\v) = k}. (58) 

A multiple-line graph T is called 1-multiple-line irreducible (1MLI) if it satisfies 
the following criterion. Decompose an arbitrary multiple-line m G A4r- Every 
remaining connected component has then at least one external line attached. It can 
be a 0-line or a [/-line. We write 

O l JfM( L ) = { r e Gh^iL) I r is imli}. (59) 

The renormalized multiple-line moment graphs are graphs that are 1LI and have 
precisely one external multiple-line, but no external vertex, 

R k (L) = {re Go L k(L) | there is m G M r with E T U \m) = k}. (60) 

The equivalence classes of graphs that are both IVI and IMLI are denoted by 

s El , E2 (L) = Gh v jE 2 (L)ng^ 2 (L). (6i) 

With the renormalized moment graphs as defined above, the lLI-susceptibilities 
are now obtained in the form 

xgi 2 = E( 2 ^) L E (62) 
l>o res El , E2 (L) 

The weights w(T) are given as a product of factors as described in the last subsection, 
with the following two exceptions. 

• The vertex coupling constants v° c are replaced by the renormalized vertex 
moments 

l>o reQ„(L) 
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• The multiple line coupling constants mjf are replaced by the renormalized 
multiple line moments 



In the series representations above, the w(T) are computed according to the rules 
of subsection 3.2. 

5 Graph construction 

In this section we describe a mechanism to generate the multiple-line graph classes 
<Se 1 ,e 2 {L), Qk{L) and TZk(L) with L > 1 internal lines. Similar as in LCEs the idea is 
to define appropriate classes of vacuum graphs (no external lines). The generation 
of the graphs of S, Q and TZ is then done by attaching external lines in various 
ways. By defining appropriate recursion relations, the order by order construction 
can then be restricted to the vacuum classes. 

Fortunately we can profit from graphs that have been generated already in stan- 
dard LCEs. This way we completely avoid a cumbersome recursive construction, but 
obtain the vacuum multiple-line graphs by operating on the LCE vacuum graphs of 
the same order. 

The graphs of S, Q and 1Z will be obtained through a sequence of simpler graph 
classes. At each step, any method of multiple-line graph construction should satisfy 
the following two conditions. 

• It should be surjective or complete. At least one graph of every equivalence 
class should be generated. 

• For every equivalence class precisely one representant should be kept. All gen- 
erated graphs of a particular class of graphs should be mutually inequivalent. 
This requires a so called non-equivalence test for every newly generated graph. 
The graph should be kept if and only if it is not equivalent to a graph already 
generated. In general it must be compared with all multiple-line graphs that 
have been generated before. If this is necessary, the non-equivalence test is 
rather time consuming. 

The algorithm we describe below circumvents an extensive comparison of mul- 
tiple-line graphs . It is defined in such a way that two graphs can be equivalent only 
if they have their origin in the same graph of the prior graph class. 

Essentially we proceed in the following steps. The starting point are the vac- 
uum diagrams of V<i{J->) as computed in [[J. Vi{V) is the set of (standard) LCE 
equivalence classes of connected, 1LI vacuum graphs with L lines. Operating on 




(64) 



L>0 



r&R v (L) 
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V2{L), we generate all mutually inequivalent multiple-line graphs with L lines that 
are connected, 1LI and have no external lines (neither 0-lines nor [/-lines), 



MV 2 (L) = Q$(L). 



(65) 



This is done in two steps. In a first step lines with the same endpoint vertices are ar- 
ranged in all possible ways into multiple lines. We obtain a subset of AiV2{L) which 
we denote by J\AVi{J^)- A multiple-line graph T G AiV2{L) belongs to AiV2(L), 
T G AiV2{L), if and only if all bare lines of T that belong to the same multiple 
line have the same endpoint vertices. The second step is to relax the additional 
constraint on MV2{L). This is achieved by resolving the vertices. This means that 
every vertex is split into several ones, and that the bare lines that were attached to 
the original vertex before now are distributed over the new vertices in all possible 
ways. As a general constraint on the resolutions, the resulting multiple-line graphs 
should be 1LI in order to stay in MV2(L). We impose an additional constraint on 
the resolutions, allowing only for so-called admissible vertex resolutions. They are 
defined in such a way that two multiple-line graphs of AiV2{L) can be equivalent 
only if they originate from the same graph of AiV2{L). This considerably simplifies 
the non-equivalence test. 

Once we have obtained AiV2(L), the final multiple-line graph classes Se!,e 2 (L), 
Qk(L) and TZk(L) are obtained by attaching external 0-lines and [/-lines in the 
appropriate way to all graphs of MV2{L). The steps in the generation of the DLCE 
graph classes are shown schematically in Fig. 6. 



Figure 6: Modules in the generation of DLCE graph classes. The root of the in- 
heritance tree are the LCE vacuum graphs of V2(L). By means of multiple-line con- 
struction (.MP2) and vertex resolution (M.V2) the DLCE vacuum graphs of AiV2{L) 
are obtained. Attaching external 0-lines and external [/-lines in well defined ways 
yields the renormalized vertices (Qfc(L)), the multiple-line moments (TZk(L)), and 
the 1MLI and 1VI graphs (S Ei ,e 2 (L)). 

In the following we describe the multiple-line graph construction in detail. Below 
we use the following notation. For any finite set A let V(A) denote the set of 
partitions of A into mutually disjoint nonempty subsets of A. For every II G V(A) 
we have 




Qk{L) 

n k {L) 




(66) 



Pen 
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5.1 V 2 (L)^MV 2 (L) 

Starting with the LCE vacuum graphs of Vi{L\ we first of all arrange internal lines 
into multiple lines. Consider an arbitrary graph T G V^iV). Simultaneously for 
every pair of neighbouring vertices v , w G £>r choose a partition II G V(C VjW ), where 
£ VtW is the set of common lines between the vertices v and w. Every / G C ViW belongs 
to a unique P G II G V(£ VtW ). For every P G II we introduce a multiple line m(P) 
such that all / G P are precisely the bare internal lines that belong to m(P). This 
procedure defines the incidence relation for all I G £r and, along with that, a 
multiple-line graph that belongs to M.Vi{p). 

In Fig. 7 we give examples for this construction. In the first one, an LCE graph 
leads to 5 mutually inequivalent DLCE gaphs, in the second one to 3 DLCE graphs 
of MV 2 (L). 

Proceeding this way for all partitions of lines for all pairs of neighboured vertices 
of T, and for all T G T > 2{L)- I we obtain the class of multiple-line graphs MV2(L). 
The non-equivalence check for every newly generated multiple-line graph T can be 
drastically restricted. It is sufficient to compare T only with those multiple-line 
graphs that have been constructed before from the same graph T G V 2 (L) as T 
results from. Other graphs of M.V 2 (L) are inequivalent to T. 




OO-OO +<|X> +<X|> 

Figure 7: Two examples of multiple-line construction V2(L) — > M.Vi{J-<) with 
L = 4. Common lines between two vertices are arranged into multiple lines in all 
possible ways. Only mutually inequivalent graphs are shown. 

5.2 MP 2 (L) -> MV 2 (L) 

The second step is to relax the constraint on MV 2 that bare lines which belong to 
the same multiple line must have the same endpoint vertices. This is achieved by 
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vertex resolutions of all graphs of MV 2 (L). To efficiently reduce the non-equivalence 
check of new graphs we allow only for so-called admissible vertex resolutions. 

Let T G AiV2(L), v E Bf any vertex of V and II G V(C V ) any partition of the 
set of lines C v entering v. We remove the vertex v and draw for every P G II a new 
vertex v (P) so that all lines / G P enter the vertex v(P) rather than v before its 
removal. This procedure is called a vertex resolution of v. The resulting graph T 
must be 1LI in order to belong to MV2{L). 

The vertex resolution of v is called admissible if in addition the (standard LCE) 
graph T is connected with T defined in the following way. For every I G C v there 
is precisely one m G Aif with (/) = m, and there is exactly one P G II with 
/ G P. Draw a vertex w{m) = w(^f(l)) for every such multiple line, a vertex v(P) 
for every P G II, and define a line /(/) with incidence relation 

$(/) = («,(%(*)), «(P)). (67) 

This defines a graph T = (Bf, £p, Ef = 0, <3>) with 

B f = {v(P)\P G 11} U {w(^f(l))\l G £ v } 

C f = {T(l)\leC v }. (68) 

A vertex resolution is called admissible if the graph T is connected. Roughly speaking 
this implies that no vertex drops off after vertex resolution. 

Fig. 8 shows two admissable and one non-admissible vertex resolutions (the aux- 
iliary graphs T are not displayed). 




Figure 8: Example of admissible ((a) and (b)) and non-admissible (c) vertex reso- 
lutions. In particular, the case of the trivial resolution (a) is admissible. 

The class A4V2(L) is now generated as follows. Consider an arbitrary graph 
r G MV 2 (L). Simultaneously for all vertices v G £>f apply an admissible vertex 
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resolution. We obtain a multiple-line graph T £ MV 2 (L). Doing this for all possible 
resolutions of all vertices and for all T e AiV2{L), we get the complete graph class 
AiV2{L). Again, due to the restriction to admissible vertex resolutions, the non- 
equivalence check of a newly generated graph T e AiV2{L) can be restricted to 
graphs that have been generated from the same graph of AiV2(L) as T was. 

In Fig. 9 we show an example for the application of admissible vertex resolu- 
tions. Here one graph of M.V2(L) leads to eight graphs of A^P 2 (L), indicating the 
proliferation of DLCE graphs. Note that the starting graph was only one of five of 
Fig. 7. 




Figure 9: Example of applying admissible vertex resolutions for AiV2(L) 
A4V2(L) with L = 4. Only mutually inequivalent graphs are shown. 



5.3 MV 2 (L) -> S Ei ,e 2 (L), Q k (L) and R k (L) 

The final step is to construct the classes of multiple-line graphs Se 1} e 2 , Qk and R k 
out of the vacuum multiple-line graphs of MV2- This step is realized by attaching 
external 0-lines and external [/-lines in different ways to every multiple-line graph 
r E MV 2 {L). 

Rk{L) is obtained by attaching k external [/-lines to just one multiple line m e 
Air, for all multiple lines m and for all T. 

Qk(L) is obtained by attaching k external 0-lines to just one vertex v e Br, for 
all vertices v and for all T. 

<Se 1 ,e 2 {L) is obtained by attaching E\ external 0-lines to the vertices of T and E2 
external [/-lines to the multiple lines of T, for all T e AiV2{L). This is done under 
the constraint that the resulting multiple-line graphs have to be 1VI and 1MLI. 

Again, the non-equivalence test of multiple-line graphs can be confined to pairs 
of graphs that have their origin in the same graph of MV2(L). 

In Fig. 10 we illustrate the attachment of external lines to a vacuum DLCE graph 
oiMV 2 {L). 
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Figure 10: Attaching external lines to a vacuum DLCE graph of AiV2{L). The 
example shows from left to right: the generation of a graph of Sq^{L), Q2(L) and 
R 2 (L), with L = 4. 

6 Applications to Spin Glasses and Neural Net- 
works 

In this section we indicate applications of DLCEs to spin glasses and partially an- 
nealed neural networks with "slow" interactions and "fast" spins Jl2|]-[|14]], jl7| . 

6.1 The SK model and the replica trick 

For simplicity we consider the Sherrington-Kirkpatrick (SK) model which is a spin 
glass model with infinite connectivity |]12|| . Following the notation of section 2, the 
fast spins are represented by the $- field $j, i = 1, . . . , N, taking values in {— 1, +1}. 
TV denotes the total number of spins. The slow interactions are mediated by the 
[/-field [/(jj), i < j = I, . . . , N. The SK model is described by the partition function 



with M some normalization that will be specified below. The effective Hamiltonian 
H of U is given by 



where \i is a positive coupling constant. Zp{U) is the partition function of a spin 

system 




(69) 




i<j=i 



(70) 




(71) 



{*i=±i} 



i<j 
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at frozen spin-spin couplings f7(ij) and temperature controlling the fast spin 
fluctuations. In (|69D it has been assumed that the equilibrium distribution of the 
slow variables is a Boltzmann distribution governed by a second temperature /J' -1 . 
This assumption is justified, if the time evolution of the Us is determined by a dissi- 
pative Langevin equation (for the precise conditions see e.g. flTj|). Such a Langevin 
equation can be actually derived for the Us from an ansatz which is motivated by 
neural networks |12|-|14 . 



One of the major quantities of interest is hxZL,. Let us first rewrite Z'o, in the 
form 



= [[Zp(U) n )l (72) 

where we have introduced J = f3fi and n = j3' / (3. The normalization has been 
chosen such that [[1]] = 1. 

The limiting case n = with J kept fixed corresponds to a quenching of the 
!7s in a spin glass, in this limit the U degrees of freedom decouple from the spin 
dynamics, while n = 1 corresponds to a pure annealing in the disordered system 
(see e.g. ||13|| ). Partial annealing then refers to a value of n with < n < 1. In 
particular the limit n — > with fixed J and /3 amounts to the replica trick in spin 
glasses |T?1. Assuming that averaging [[■]] and limes n — ► commute, we have 

[[In Z p {U)]) = lim^k (T3) 

The left hand side of (1731) is calculated by means of the right hand side. 

So far n has been a positive real number. For integer n = 1, 2, . . ., Zp(U) n is the 
partition function of an n-times replicated system. Let us first consider this case. 
We rewrite 

n N 

ZpiPT = £ ex P ((3j2 E Vto)*™*?)' ( 74 ) 

a=1 i<J=l 

with a = 1, . . . , n labelling the replicated spin variables, so that 

poo N 

Z 'np = / U dU ^)- E exp(-S(U,QM)), 

n AT 1 
a=l i<j'=l i<j 

Linear terms in $ and {7 may be included according to 

n N N 

s lm = -^EE $ S a)+c E tfw) ( 76 ) 

a=l i=l i<j=l 
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with constant external fields h and c. Apparently Z' n * has the form of models to 
which DLCE applies, with a hopping term 

n N 

s hop (u,$ (a) ) = -/^EE u«j)*P*P> ( 77 ) 

a=l i<j=l 

a single link action 

SWm) = cU {h3) + l -JNUl jy (78) 

and a single site action 

n 

S°(^) = -hY,^- ( 79 ) 

a=l 

Application of DLCE to In Zp, and derived quantities yields their series expansions 
in (3. 

So far, merely for simplicity, we have considered the SK model which is a spin 
glass model with infinite and uniform connectivity, in which all spins are coupled 
with equal variance [[f^]] = 1/(JN). For integer n — 1, 2 . . ., this system becomes 
exactly solvable in the thermodynamic limit corresponding to the limit N — > oo, so 
that the leading order of a large- N expansion becomes exact, cf. e.g. [|I7[]. In DLCE 
this amounts to a resummation of appropriate tree graphs. 

In cases of finite or non-uniform connectivity it is no longer possible to obtain 
solutions in a closed form. Examples are multilayered neural networks with couplings 
between adjacent layers, and spin glasses with short range interactions. Models of 
this kind are suited for an application of DLCEs and will be studied in future work. 

Finally we come to the limiting case of n = 0. Usually it is obtained by extrapo- 
lating the results obtained for n — 1, 2, . . . towards n — 0. For instance, [[In Zp(U)]] 

is obtained from [[Z ( g(?7) n ]] for n — 1, 2, This is a rather subtle point because the 

extrapolation is not unique without further assumptions. In the next subsection we 
show how it is possible to completely avoid both the extrapolation and the replica 
trick and to directly compute the quantity of interest by means of dynamical linked 
cluster expansions. The task is to appropriately identify the classes of contributing 
graphs. 



6.2 Avoiding the replica trick 

First we adapt the notation to section 2 to include more general cases. A denotes 
the support of the spins, that is the set of lattice sites, with V = |A | denoting their 
total number. Ai C Ai are the pairs of sites whose spins interact. In accordance 
with (|72|), we write for the normalized link-average of a function f(U) 

[[/([/)]] = Jw f(U) (80) 
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with 

VU = Y[dn(U(l)), 

/oo 
dfi(U) = 1. (81) 
-oo 

It is convenient to introduce the single link expectation values 

< 9(U) >i = j dn{U) g(U) (82) 

and the generating function W /1 (7) by 

expW\l) = < exp (IU) >i . (83) 

The way in which the replica trick can be avoided is examplified for the free 
energy density W sp /V of the spin system averaged over the link couplings. The 
partition function of the spin system for a given distribution of the link interactions 
U(x, y) is given by 

ex P W sp (U) = N sp J V$ exp(-S sp ($,U)), (84) 
where W sp (0) = and 

S sp ($,U) = -\ v(x,y)$(x)$(y)U(x iy ), 

x,y€A 

p$ = Y[d$(x)-exp(-S°mx))). (85) 

x£A 

Without loss of generality we identify the support of the interaction v(x, y) = v(y, x) 
with the set A x of lattice sites where VU is supported, 



A x = {I = (x, y) e A x A | v(x, y) ^ 0}. (86) 

For simplicity we assume v(x,y) to be of the form fl46|) , so that K is a measure of 
the strength of the interactions v(x,y). 

The free energy density of the spin system allows for a series expansion in the 
standard LCE sense, with the link field U (I) playing the role of a "background field" , 

^w sp (u) = J2( 2K ) L E ^ sp (r,f/)- (87) 

l>o reg s p (L) 
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Here Q S ^{L) (with E = 0) denotes the set of equivalence classes of connected LCE 
graphs with E external lines and L internal lines. The spin-weights w sp (T, U) are of 
the form 



w 



sp 



(T,U) = R sp {T) J2 YlU(l) m ®. (88) 

The sum is taken over all non- vanishing lattice embeddings of the graph T. It runs 
over all maps of internal lines of the graph T to pairs of lattice sites of Ai that are 
consistent with the graph topology in the sense discussed in section 3. For every 
I £ Ai, m(l) denotes the number of lines of F that are mapped onto the link I by 
the embedding. All other factors that contribute to the weight are collected in the 
prefactor R sp (T), including the inverse topological symmetry number of T. 

Next we want to express [[W sp ({/)]] as a series in K by means of DLCE. Toward 
this end we set f(U) = W sp {U) and insert the series (|87| ) with (|88|) into (|80|). At this 
stage we are not concerned with question of (uniform or dominated) convergence 
and obtain 



[[^ W s P 



(U)}} 



L>0 



E 

ree sp (L) 



VU w sp (T, U) 



RSP ( T ) Yl II < u ( l ) m{l) >i 



(89) 



L>0 



The next step is to express the full single link expectation values in terms of the 
connected ones. They are related by 



< U m >i 



e n 

nev(m) Pen 



< 



17' J 



> 



i- 



(90) 



where V(m) denotes the set of all partitions of m = {1, . . . , m} into non-empty, 
mutually disjoint subsets of m. \P\ is the number of elements of the set P. The 
relation ( j90"D is equivalent to the partition of all lines of T that are mapped to the 
same lattice link into multiple lines, with every multiple line contributing a factor 



Using @g), (0) we rewrite as 



d\ p \w\i) 



dP 1 



m 



1=0 



lc 

pi- 



rn 



d w sp (u))) = J2my 

L>0 



ree ( f(L) 




e ri 1 

Pen / \n^Ai leAi 

' (92) 

The last summation in fl92l ) is over all maps Cr Ai of the lines of T to the 
lattice links of Ai subject to the constraint that all lines that belong to the same 
multiple-line corresponding to some P £ II are mapped onto the same lattice link. 
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Finally we rewrite (|9"2"| ) as a sum over multiple-line graphs . To this end, we 
first observe that for every T G Qq P (L), every partition II G V(Cr) of the lines 
of r into multiple-lines generates a multiple-line graph A = (T, II) in the obvious 
way. Let us denote by Q 0t0 (L) the subset of multiple-line graphs of Q 0i0 (L) that 
stay connected after decomposition of all multiple lines. (These are the multiple- 
-line graphs which stay connected in the usual graph theoretical sense, when the 
dashed lines are omitted.) For every A G Q 0t0 (L) there is a unique T(A) G Qq P (L) 
and at least one II G V(Cr(A)) such that (r(A),IT) = A. Let ba be the (uniquely 
determined) number of partitions II G "P(£r(A)) with (r(A),II) = A, and n(A) 



such an arbitrary partition. Eq. ( |92| ) then becomes 



4 W sp (U))} = J2( 2K ) L E ™A^(r(A)) n mfa E_ 1 

l>o Aeg 00 (L) \Pen(A) / \n(A)->Ai 

(93) 

The last bracket of (|93D is the lattice embedding factor of the multiple-line graph 
A. The second bracket from the right does not depend on the choice of n(A) and is 
the product of the multiple-line coupling constants as defined in section 3. Finally, 
?7.A-R sp (r(A)) is precisely the remaining part of the weight of A that was described in 
detail in section 3, endowed with the correct inverse topological symmetry number 
of the multiple-line graph A (because of the factor ri/\). 

In summary, we obtain the series expansion of the link-averaged free energy 
density in terms of DLCE graphs, 

4 w »w]] = E( 2K ) L E ( 94 ) 

L>0 Aeg , (L) 

The weight w(A) of a multiple-line graph A is defined and computed according to 
the rules given in section 3. 

Eq. ( p4|) is the series representation of the link-averaged free energy density of 
the spin system, i.e. the free energy density of the n = replica system, in terms of 
DLCE graphs. It looks much like the series representation of the 1-replica system, 
which is given by 

^W 1 . repl = i \n[[ex P W sp (U)]] = E^ E W ( A ) ( 95 ) 

L>0 Aeg Q , (L) 

according to the discussion of section 2. We recall that Q fi(L) is the set of DLCE 
vacuum graphs with L bare lines that are connected in the generalized DLCE sense. 
Comparing (0) and (|95|) , the transition from n = 1 to n = replicas is achieved by 
keeping only the subset Q ^(L) C £ (L) of multiple-line graphs that are connected 
in the original (LCE) sense. 

We expect that the series (0) are convergent for a large class of interactions 
S l (U) and v(x,y) if the coupling constant K is sufficiently small. For special in- 
teractions most of the multiple-line graphs yield vanishing contributions so that we 
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can further restrict the sum to a subset of Go,o (L). An example is given by the mean 
field type of interaction of the SK model where all pairs of sites are coupled with 
the same strength. That is, in the notation of this subsection, 

v(x,y) = K(l-S x , y ), 

Ai = Ai is the set of all pairs of sites, (96) 

sHu) = v-u 2 . 

v ' 2 

In the thermodynamic limit V — > oo only those multiple-line graphs survive that 
consist exclusively of 2-lines, and that do not contain any loop of 2-lines. The 
resulting tree structure provides the possibility of summing the series (|9"3|) and of 
analytically continuing the result to the spin glass phase of large coupling K. Work 
in this direction is in progress. 



7 Applications to the SU(2) Higgs model 

In this section we discuss applications of DLCEs to study the Higgs transition within 
a variational estimate for the gauged SU(2) Higgs model. The estimate serves as an 
effective description of the electroweak standard model. 

We consider a 4-dimensional hypercubic finite temperature lattice Ao of size 
Lq x V3, with Lq = T~ l the inverse temperature in lattice units and V3 the spatial 
volume. The lattice links are the set of nearest neighbour lattice sites, given by 

A x = {(x; fi) I x G A , n = 0, . . . , 3}. (97) 

The gauge field U(x; /1) is an SU(2) valued field living on the lattice links Ai. It is 
convenient to parametrize such an SU(2) matrix by 

U = 4> (U)l 2 + ia-$(U) 

= U + ia-U (98) 

with U ~ (Uq, U) G 6*3. The Higgs field <3> lives on the lattice sites A . Its values are 
real multiples of SU(2), 

= <f>o{&(x))l 2 + ia-$($(x)) 

= (j> (x)+ia-$(x), (99) 

with ~ (<p (x) , <j)(x)) G R 4 . We say that $ is cU{2) valued. 
As a special case of Eq. (|I|) we consider the partition function 

Zvei = tvUV® exp(-S VE i(U,$)), (100) 
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in which 

3 

VU = II HdVniUfcfi)), (101) 

with dfXniU) the normalized Haar measure on SU(2), and 

X?$ = Y[ exp(-5°($(x))) dv(<S>(x)), (102) 



where 



dz/($) = rf 



i. 



= itr($t$) + A Qtr($t$)_l^ 
= 2 + A(0 2 -l) 2 . 



(103) 



Finally, 



5y B i(C/,$) = -^{4O i n fe 0o(t/(^;O)) + ^4C cu6e 0o(t/(^;^)) 

3 1 

+ ^(2«)-tr (^(^[/(x; + £)) (104) 

+e0o($(x)) 

depending on variational parameters Clink, (cube and £ G R, and on the hopping 
parameter k. The first two terms are non-gauge invariant substitutes for the Wilson 
plaquette term of the gauge part of the SU(2) Higgs model, cf. Eq. ( |106| ) below. We 
distinguish timelike from spacelike directions. Timelike links are denoted as (x; 0) 
and spacelike as (x; /i), \i = 1, 2, 3. The third term is a hopping term coupling $ and 
U fields only in spatial directions. The last term depends on the third variational 
parameter £ associated with the Higgs field $, while the remaining ultralocal action 
for $ has been absorbed in the measure. 

In the next subsection we first indicate how we arrive at the action ( |104| ) within 
a variational estimate for the free energy of the SU(2) Higgs model. We derive 
equations for the optimal choice of variational parameters and an equation -based 
on stability arguments- to determine the critical hopping parameter K crit at which 
the Higgs phase transition sets in. These equations will depend on expectation 
values of n-point correlations that are most suitably evaluated with a DLCE in 
three dimensions. 

The generic graphical expansion exposed in sections 3-5 now has to be adapted to 
account for the internal symmetry. This is somewhat involved so that we restrict our 
discussion to some remarks in section 7.2, but postpone any details to a forthcoming 



paper [fT9 |. In section 7.3 we present some results for the critical line in comparision 



with other variational estimates and Monte Carlo simulations. 
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7.1 Variational estimates for the free energy of the SU(2) 
Higgs model 

By means of DLCEs we want to obtain an analytic estimate for the critical line in the 
SU(2) Higgs model that serves as an effective description of the Higgs transition in 
the electroweak standard model. Because of the rather different methodical ansatz 
in this approach, for a common choice of parameters this estimate will provide 
an independent check of numerical results based on Monte Carlo simulations. In 
addition, larger scalar field couplings are available in our approach and larger values 
for L , the number of time slices. (In our calculations L = 4 or L = oo, while it 



was chosen as 2 or 3 in |15|.) Larger Lqs are required for a finite size scaling control. 



To fix the notation we introduce the action of the SU(2) Higgs model as 



Z = J VUV<$> exp (-£(£/,$)), 
S(U,$) = S w (U) + S hop , (105) 

3 j 

S hop = -^{^(2 K )-tr($(x)tf/(x ;/ ,)$(x + /i))}. 



The Higgs self- interaction has been absorbed in the measure P<3> of Eq. (102), Sh op 
is the hopping parameter term in 4 dimensions. Sw is the SU(2) gauge invariant 
Wilson action 

S W {U) = /Z P f 1 ~ ^^U(x]^)U(x + fi]u)U(x + u]fi)- 1 U(x]v)- 1> \ 

xeAo n<v=0 ^ ' 

(106) 

with gauge coupling (3 = 4/ ' g 2 and U G SU(2) as in (|104|) . 

A plaquette term involves a 4-link interaction. In contrast to the hopping param- 
eter term it does not allow for a direct treatment with DLCEs. Therefore we replace 
Z, the full partition function of the SU(2) Higgs model, by a partition function Zve 
that is related to Z by an inequality of the form 

exp (-Vf) = Z > Z VE exp (< - (S - S VE (Q) >ve) = exp (-Vf(Q). (107) 

Eq. ( |107| ) follows from the convexity of the exponential function and holds indepen- 
dently of the specific choice of Z and Zve, if the measure is positive definite and 
normalized. V is the 4-dimensional volume, / denotes the true physical free energy 
density defined via Z, f the trial free energy density defined via the second equality 
in ( 107 ). S refers to the original action in Z, in our case it is the SU(2) Higgs action, 



and Sve to an ansatz for the action in Zve depending on a generic set of variational 
parameters (. These parameters should be optimized so that the difference of the 
trial and the physical free energy density f(Q — f > becomes minimal for £ = £ 
at which / takes its minimum. 
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The most naive ansatz for Sve is a mean field ansatz in the spirit of a molecular 
field approximation, leading to a factorization of Z VE0 according to 

z VE0 (c,o = zf^(0-z^(0. (108) 

The partition function depends on two variational parameters ( and £ and factor- 
izes in a product over single link (Zi ink ) and single site (Z site ) partition functions, 
V3 denotes the three-dimensional volume. In this ansatz there is no space for imple- 
menting an asymmetry between temporal and spatial directions. Thus the results 
will be temperature independent by construction. 

One can think of a variety of improvements of the molecular field ansatz, argu- 



ments against or in favour of the various versions will be given in ||19|| . Here we only 
consider the ansatz which is most plausible from a physical point of view and leads 
to the best agreement with Monte Carlo results for comparable sets of parameters. 
The ansatz treats the spacelike degrees of freedom of the hopping term in 3 dimen- 
sions beyond the mean field level (that is with DLCE), but the timelike degrees of 
freedom within a mean field approach, so that the partition function factorizes over 
the spacelike hyperplanes. The reason is that the spatial hopping term is supposed 
to contain the nonperturbative properties of the full model that drive the Higgs 
phase transition. While a mean field approach for all variables will be too rough to 
produce high quality results for n C rit, it appears more reasonable for timelike vari- 
ables, although a finite temperature effect on the Higgs transition gets lost this way. 
(We should remark that from results in a scalar $ 4 theory in four dimensions the 
finite temperature effect on n crit is expected to be anyway quite small j|, H.} 



Now it is easily checked that the choice of Svei in Eq. ( |104| ) leads to a factor- 
ization of Zvei according to 

ZvEl((linki (cube, = Z cu b e ((cube, £) ° ' ^Unk(Clink) L ° 3 ■ (109) 

The partition function Z cube of the spacelike degrees of freedom reads 

Z cube = I J] (d^x))e- s °^]Jdfi H (U(x;^)\ • exp (-S cu6e ([7, $)), 

(110) 



(3) 

with Aq the 3-dimensional lattice and 

3 



3 ^ 

+ ^(2/ t )-tr($(x) t f/(x;/i)<l>(x + /2))}, (111) 



2 

At=l 



while Zi in k is an ultralocal one-link integral 

J dfi H (U) exp(-S Hnk (U)) (112) 



Hink 
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with 



ShnkiU) = 4:( Hnk (j) Q (U). 



(113) 



Expectation values < O > cu be, < O >n n k, < O >vei of observables O refer to 
Z cu b e , Zu n k and Zvei, respectively. Minimization of the trial free energy density 
f (Clink, (cute, m terms of these expecation values leads to three equations. The 

first one (df/d£j = is solved by £ = in the symmetric phase. The remaining 
two equations are given by 



cube 



4 dC 



cube 



(3{< (f) {U) > cuhe f + f (< MU) >linkf < MU) >cube -Ccube 







< MU) > 



cube 



(114) 



cube 



and 



3(3 



(< <f> Q (U) > cube ) 2 < <j> Q (U) >l ink -Clink 







dC 



< 4> {u) > Hnk 



0. (115) 



link 



Here we have used the shorthand notation 



< <f) (U) > cube = < 4> (U(x; 1)) > 



cube, 



^IL = \ < tr U{x; l)U{x + T; 2)U{x + 2; l)- x ?7(x; 2)" 1 > cu6e (H6) 

- {< MU{X;1)) > cu be) 4 , 

with x G Aq 3 \ In (|116|) we have used the lattice symmetries. Eq.s ( |114j) , ( |115|) 
should be solved for Clink, Ccube, and £ as series in n. The stability condition for the 
symmetric minimum at £ = is given as 



3(3 



+ 12 



d 2 W}' 2 



cube 







(3(< (f> (U) >cube) + 7T < MU) >cube (< M U ) >link) ~ C 



cube 



2 < (j)o(U) >cube 



+ 
< 



?=0 





(117) 



(2k) 2 < <f> (U) > Unk — < O ($) > cube 



- 1 



^- < O ($) > cu be 

del 



with 



< 00 ($) >c«be = < 0o($O)) > 



cube, 



;ii8) 



a; G Aq 3 ^. For an equality sign, (|117|) determines K cr j t order by order in the expansion. 
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Note the two type of expectation values < • >u n k, < • >cube entering Eq.s 
( |115|) , ( |117|) . While the single link expectation values < • >u n k can be evaluated 
exactly, the < ■ > C Mfe e s must be approximated. The derivatives of W c ^ he w.r.t. Ccube 
or twice w.r.t. £ induce up to 5-point functions in U and 6-point functions in four Us 
and two $s. Because of the bad signal/noise ratio Monte Carlo calculations of such 
connected correlations would not be feasible. Therefore this is the place, where we 
use a DLCE to evaluate expectation values of the type < [/•••?/<&•••$ > cu be (all 
internal and configuration space indices suppressed) order by order in k. Clearly, 
also in this scheme the price is high. For example, the number of connected DLCE 
graphs contributing to order k 4 to a 4-point function of 2 Us and 2 $s is about 
100. Not all correlations in ( 114| ),( |TT5D ,( |117|) appear in the form of susceptibilities. 
In a product of 4 Us, for instance, the configuration space indices are fixed to the 
boundary of a plaquette. Such features must be noticed for calculating the graphical 
weights. 

These remarks may indicate the complexity of the actual evaluation of Eq.s (114), 
( 115 ), ( 11 7]) . Here we omit any further details on the list of contributing graphs 
and postpone them to [PL9|| . Before we quote some results, we comment on further 
subtleties, when DLCEs are adapted to internal symmetries of the hopping term. 



7.2 DLCEs for internal symmetries 

In this section we focus on the (gauge) group SU(2), since we are interested in 
applications to the electroweak phase transition, but it should be obvious, how 
one could proceed along the same lines for a group SU(N) with N > 2. In the 
development of the multiple-line graph theory in sections 3-5 we had assumed that 
(f) G R, U G R to simplify the notation. In Eqs. (|9|) and (fH) we had $ G cU(2), 
U G SU(2). Let us write the hopping part of the action in the symmetrized form 



Shop = -~ v(x,y)ltr<S>(xyU(x,y)$(y), (119) 
with hopping parameter 



2 

x,y£Ao 



3 



v(x, y) =2k^ ($y,x+n + 5 y,x-Ji) ( 12 0) 

At=0 

and 

U(x, x + fi,) = U(x;fi) 

U(x + fi,x) = U(x;fi)^ 1 , fi = 0, . . . , 3, (121) 
U(x, y) =0 otherwise 



Using the parametrizations (|98|) and (|99]) we have 
I 

-tr$(x)W(:r,y)$(y) = <f>(x) ■ <j>{y) U (x, y) 
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+Mv) <l>(x) ■ U(x, y) - o (ar) <j>{y) ■ U{x, y) (122) 
+4>{x) x $(y) -U(x,y). 

Note that in the parametrization in terms of (<po, <j))s and (Uq, U)s the S° part 
of the action is 0(4) invariant, but Svei is only 0(3) invariant for ( ^ or £ ^ 0. 
Thus we have to distinguish between singlet (0) and triplet (s = 1,2,3) indices in 
internal space. 

Recall the generating equations @ of the graphical expansion of a DLCE in 
which < Q(x)U(x,y)Q(y) > had been expressed in terms of connected n-point 
functions, n — 1, 2, 3. Using the 0(3) symmetry it follows from Eq. ( 122 ) that each 
derivative d/dl(x,y) now has to be replaced by the differential operator 

d 3 9 

Lafe(X ' y) = **aZbfri) + ^ (5a ' A ° " + ^ dL^y)- (123) 

Latin indices run from to 3, Greek indices from 1 to 3. Furthermore we have 
introduced the totally antisymmetric e-symbol which is zero if any of the three 
indices is zero and €123 = 1. This implies for instance for the one nonvanishing term 
at v = in Eq. @ a substitution according to 

W H{x) W H(y) W I{Xty) — ► W„a (x) W Hb(y) L ab (x,y)W. (124) 

According to the graphical rules an ra-line was generated by n iterated derivatives 
w.r.t. Ij,j G 0, 1, 2, 3, evaluated at v — 0. These derivatives now have to be replaced 
by a multiple application of L yielding an n-line 




3 



where the 1-link generating function W l (I) is defined by 



(125) 

7=0 



P ^ /^(C/) exp(4C CMfee o (^)) 1 ' 

For cii = hi — 1 for all i G l,...,n, n applications of L a b to IU 1 reduce to 
5 a bd n W 1 1 dl n 1 7=o, but otherwise we get a complicated mixing of singlet and triplet 
terms. 

To evaluate products of Ls systematically, it is again quite convenient to manip- 
ulate graphical rather than analytic expressions. We represent one application of 
L ao on W 1 as shown in Fig. flTTp. Multiple applications of L a b require an iteration 
of Fig. (0). 

Twiddle lines are bare 0-lines with a singlet index a = 0, crossed lines are bare 
0-lines with a triplet index a = 1,2, 3. The arrow in the last term stands for the 
e a fe 7 symbol in L ab , its orientation reflects the order of indices airy. 
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Figure 11: Graphical representation of L^W 1 . It provides the generating equation 
for the graphical computation of the internal symmetry factor of a multiple-line. It 
replaces Fig. 1(b) for / = H = 0, v = and n = 1. For further explanations see the 
text. 



This way one achieves a further graphical "decomposition" of an n-line into 
terms merely consisting of bare lines, from which the internal symmetry factors can 
be easily read off. For example a 2-line is further decomposed into 11 terms, a 3-line 
into 29 and a 4-line into 142. This may indicate that the number of decomposed 
contributions rapidly increases with the number n of bare lines collected to an n-line. 
In our actual calculations up to order k 4 the maximal number of internal lines was 
4 and the highest multiple line was a 6-line with 4 internal and 2 external U-lines. 
Some of the results will be reported in the next section. 



7.3 Some results 

In Table 1 we compare results for K crit obtained for 2 values of f3 and A and various 
lattice volumes within three methods. The first line refers to a naive molecular field 
approximation with a partition function Zveo factorizing according to Eq. ( |108| ). 
Because of the complete factorization in space and time directions the ansatz keeps 
only a trivial volume dependence and is temperature independent by construction. 
The second line refers to Monte Carlo simulations of the SU(2) Higgs model in 



four dimensions [15| for two lattice extensions Lq = 2 and 3. 

The third line shows results of the DLCEs with the variational ansatz (VE1) 
presented in section 7.2. Again this ansatz is independent of Lq by construction, so 
that we get an identical result for Lq = oo. The values for K crit have been obtained 
by a linear regression in 1/R 2 , R denoting the order in k of the DLCE. For the 
regression we used the coefficients of k 1 , k 2 , k 3 and k 4 . The value of k ct h at order k, 4 
was given by 0.1284 for A = 5 • 10 -4 and (3 = 8.0. Although the fourth order seems 
to be low for a series expansion, the agreement with the Monte Carlo results is not 
too surprising, if we remind the complexity of the expansion, which is manifest in 
the number of graphs contributing at fourth order to the various expectation values. 

Furthermore, for otherwise fixed parameters, we see a decrease of k ct h between 
the Monte Carlo and the DLCE results for increasing Lq. This is in qualitative 
agreement with the expectation. 
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Method 


P 


A 


lattice 




Mean Field 


8.0 


5.0 


• io- 


-4 


oo 


4 




0.12973 


Monte Carlo 


8.0 


5.0 


• io- 


-4 


2 • 


32 • 32 


256 


0.12887(1) 


VEl & DLCE 


8.0 


5.0 


• 10" 


-4 


4- 


oo 3 




0.1282(1) 


Mean Field 


8.15 


5.1 


• io- 


-4 


oo 


4 




0.12964 


Monte Carlo 


8.15 


5.1 


• io- 


-4 


3- 


48 ■ 48 


384 


0.12852(2) 


VEl & DLCE 


8.15 


5.1 


• 10" 


-4 


4- 


oo 3 




0.1281(1) 



Table 1: Results for K cr u for 2 values of (3 and A, various lattice extensions, and 3 
methods. The Monte Carlo results are taken from [15], the mean field and DLCE 
results, extrapolated to infinite order, from this paper. (The error in the DLCE 
results refers to the extrapolation of n C rit-) 



In Table 2 we compare results for f« cr #, to order k a (first column) and extrapolated 
to infinite order (with the same regression as explained above) (second column) for 
DLCEs with different variational estimates. The couplings are fixed, and Lq = 4. 
VEl refers to the first variational ansatz, also exposed in Table 1. The ansatze for 



the variational estimates VE2 and VE3 will be explained in detail in [19]. Here 
we only characterize them by their different factorization properties of the partition 
function. 



Ansatz 


0(/? 4 ) 


extrapolated 


VEl 
VE2 
VE3 


0.1284 

0.12958 

0.12948 


0.1282(1) 

0.12910(6) 

0.12934(7) 



Table 2: Results for K cr u, obtained in a DLCE to order k 4 (second column) and 
extrapolated to infinite order (third column) at (3 = 8.0, A = 5 • 10 -4 , L = 4. The 
ansatze VEl, VE2, VE3 refer to 3 variational estimates for the free energy of an 
SU(2) Higgs model. Their gap equations are solved with DLCE. 

The partition function ZyE2 factorizes acording to 

ZyE2 (C(0, string) i C(s,link) > £) ^(0, string) (C(0, string) > £) ' ^(s,link) (C(s,Unk)r L0V3 , (127) 

i.e., a product of a partition function Z(Q >str ing) that couples fields with hopping 
parameter terms along a string in time direction, and Zr s ,ivnk)i a partition function 
for a single spatial link. While expectation values w.r.t. Zunnfy can be evaluated 
exactly (within the numerical accuracy), expectation values w.r.t. ZcQ jString \ are 
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preferably evaluated with DLCE. This way they do depend on L , in principle they 
are sensitive to detect a temperature effect on n crit . So far the couplings along a 
string are only kept in time direction, while the spacelike directions are incorporated 
on a mean field level in Z^^unk)- 

The third variational ansatz VE3 therefore is chosen as a symmetrized version 
of VE2 in that the string now can extend in any of the 4 directions. Again we 
treat space-and timelike degrees of freedom differently, because we are interested 
in systems at finite temperature with temperature dependent coupling K cr u- The 
number of variational parameters is larger, each of it depending on two indices telling 
us whether the parameter belongs to a single link action or to a string, and whether 
the string extends in temporal or spacelike directions. 

For both versions VE2 and VE3 the actual L dependence (L = 4, oo) lies 
outside the systematic error, if we expand the series to order k, 4 . For Lq = 4 
the fourth order is the minimal order in k at which a finite temperature effect is 
visible in principle. (For the graphs it should be possible to wind at least once 
around the torus in O-direction.) At 0(/t 4 ) the part of graphs for which the lattice 
embedding is sensitive to L is about 1% of the total number of graphs. Such a 
small "signal/noise" ratio produces an effect in the 7th or 8th digit, outside the 
error of the O(10~ 4 ) because of the truncation of the DLCE series. By increasing 
the order of the expansion, finite-temperature effects will become more pronouced. 
For higher orders an algorithmic implementation of the proliferating graphs becomes 
unavoidable. 

8 Summary and Conclusions 

In this paper we have introduced a new expansion scheme for 3-point interactions 
or, more precisely, for point-link-point interactions. This scheme generalizes linked 
cluster expansions for 2-point interactions by including hopping parameter terms 
endowed with their own dynamics. In chapters 3-5 we have developed a multiple- 
line graph theory with an additional new type of multiple-line connectivity. We 
have introduced appropriate equivalence classes of graphs and discussed the issue of 
renormalization. The main building blocks for an algorithmic generation of graphs 
have been constructed. Because of the fast proliferation of graphs already at low 
orders in the expansion, a computer aided implementation becomes unavoidable, if 
one is interested in higher orders of the expansion than we have computed so far. 

In chapter 6 we have indicated promising applications to spin glass systems. 
In particular we have identified the DLCE graphs contributing to the link-average 
[[mZg([/)]] of the spin free energy \n.Zp{U). In the past such quantities often were 
only accessible by means of the replica trick. 

In section 7 we have outlined encouraging results for the transition line of the 
electroweak phase transition within the SU(2)-Higgs model. The results have been 
obtained by DLCEs applied to gap equations that follow from variational estimates 
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of the free energy. They are in good agreement with corresponding high precision 
Monte Carlo results. As we have seen, the main complications of DLCEs for SU(N) 
Higgs systems come from the internal symmetry structure of the hopping term. In a 
forthcoming paper we will demonstrate how one has to refine the graphical represen- 
tation to calculate internal symmetry factors merely within a graphical expansion. 

In conclusion, DLCEs are involved, but practicable, at least with computer aided 
generation of graphs. They provide an analytical tool to study systems in situations 
in which it has been impossible so far. 
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